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We derive the momentum space dynamic equations and state functions for one dimensional quan- 
(<-^ I tum walks by using linear systems and Lie group theory. The momentum space provides an analytic 

^— ^ , capability similar to that contributed by the z transform in discrete systems theory. The state 

fvj ' functions at each time step are expressed as a simple sum of three Chebyshev polynomials. The 

functions provide an analytic expression for the development of the walks with time. 
>»■ 
<^ 

S: 

'nJ- I I. INTRODUCTION 

(N. 

*vq ■ The study of quantum walks has received considerable attention since the introductory papers on the subject, such 

^ ' ^s [l|,Q| and references therein. In this paper, we develop an analytic approach to study the properties of these walks 
t — . based on a momentum space representation. 

0^ ■ This paper is structured such that in Section 2 of the paper the momentum space dynamic equations for one dimen- 

^~^ , sional quantum walks are derived via the Z transform of the position space dynamic equations and its representation 
^^ ■ of the discrete Fourier transform when Z lies on the unit circle. An exponential form of of the momentum space time 
vQ ' operator is derived in section 3 by using the group theory of SU{2) and a matrix inner product space. The exponential 
^^ . form allows a simple analytic calculation of the time evolution operator for arbitrary time intervals. This is used in 
■^ I Section 4 to obtain analytic expressions for the momentum space wave functions of quantum walks at arbitrary times. 
These wave functions are expressed quite simply in terms of Chebyshev Polynomials of the second kind. Some plots 
of the momentum space probability densities for different parameter values and times are provided in section 5. The 
-H , conclusions are summarised in Section 6. 

Q^. II. MOMENTUM SPACE DYNAMIC EQUATIONS 



Oh! 



^ 



For a given -0(0,0) we consider the evolution of a quantum state ip{t,x) G C^ for discrete times i > on a line 
X ^ Z. The dynamics of the state then evolve according to the difference equations, 

ipo{t> x) = e*"[a0o(i -l,x-l)+bi;i{t-l,x~ 1)], 
0i(i,2;) = e*"[-6*V'o(i - l,.x + 1) + a*^tjji{t -l,x+ 1)], (1) 

where |ap + |6p — 1 and a e R. 

Taking two-dimensional Z transforms of these equations yields 

i^o{zi,Z2) = e'" Zi^ Z2^[a^o{zi, Z2) + b'ipi{zi,Z2) 
4'i{zi,Z2) = e'^^z^^ Z2'^[-b*i>o{zi, Z2) + a*Tpi{zi,Z2). (2) 



Thus the transfer matrix for the system is 

B{zi,Z2) 



azj bz2 
—b*Z2 a*Z2 



(3) 
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therefore, for any iteration (time) index n, the quantum walk state ^(n, x) has transform x '^ z 

*(n,x)<->e''""C"(z)*(0,0), 
where C(z) is the matrix polynomial 



C(z) = 



—b * z a* z 



(4) 



(5) 



It should be noted that C is paraunitary, that is C ^{z) — C-^(l/z). In particular this impHes that C(z) is unitary 
on \z\ — 1. Further we note that detC{e'^P^) = 1 and hence the matrix 



Sip) = Cie^P) 
is unimodular. The Fourier transform x ^^ p is 

^(n,a;)^e™"S'"(p)^(0,0). 



(6) 



(7) 



Thus by choosing Planck's constant h — 1, the momentum space representation of the quantum walk state vector 
(j){n,p) evolves as 



Kp)=e™"^"(p)0(O,p), 



where 



(/.(o,p) = v(o,o) 



^o(0,0) 
^i(0,0) 



(8) 



(9) 



Thus the time evolution operator in the momentum space is a 2 x 2 matrix polynomial. Hence, the momentum 
space equations are much more amenable to analysis than those in position space. 



III. EXPONENTIATION OF THE TIME EVOLUTION OPERATOR 



The unimodular matrix S{p) can be written in exponential form as 

S{p) = Exp{ie{p)-^{p).-^) 
where 6 and 7? are real functions of p and the matrix vector if has Pauli matrix components J3| 



CTi 



1 

1 



(10) 



and 



0-2 



0-3 



-i 
i 



1 
-1 



The inner product 



iA,B)^-TriAB) 



(11) 



defined on the vector space of 2 x 2 unitary matrices gives an inner product space. The set of matrices {/, ci , 172, 0-3}, 
provide an ortho-normal basis for this space. 

The coefficients of the matrices can be evaluated by taking the inner product of both sides of (jlO|l 



ia,,S{p)) = ia,,Exp{ie{p)7^{p).^) 



with each of the matrices ai. In doing this we note that a generaUsed de-Moivre principle gives 

Exp{i6~cf .~^) = Icos{6) + ilf.lfsin{0), 
where the p dpendence has been suppressed for simplicity. Hence, 

(/, Exp{ie-^ .-^)) = cos{9) (12) 

and 

{(jj,Exp{ie~(^.~^)) = icjsin{e). (13) 

The equivalent coefHcients for S{p) can be obtained by defining 

a = cos{(3)e~''^ , 

b = sin{f3)e-'^. (14) 

Substituting in IJHJ gives 



S{p) 






(15) 



These expressions can be simplified by setting p' = p + j and p" = p + 5. Using de Moivre's principle once again 
we obtain the transition matrix coefficients 

{I,S{p))=cos{(3)cos{p'), 
{ai,S{p)) — —isin{p)sin{p"), 
{a2,S{p)) = isin{p)cos{p"), 
{(73, S{p)) = -icos{/3)sin{p'). (16) 

Comparing coefficients in equations Ijl2|l and Ijl3|l with those of Ijl6|l we obtain 

cos{0) = cos(f3)cos{p'), 

Cisin{9) = —sin((3)sin{p"), 
C2sin{9) = sin{P)cos{p"), 
C3sin{9) = ~cos{(3)sin{p'). (17) 

IV. MOMENTUM SPACE STATE FUNCTIONS 

A dynamic equation for momentum space state functions was given in (jHJ . The exponentiation of the operator in 
(fnUl enables us to write the powers of the evolution operator as 

S^{p) — Exp(in6~c? .1^) ~ Icos{n6) + i~c? .~a^ sin{n9) . (18) 

The trigonometric expressions in the above equation can be expressed in terms of the Chebyshev polynomials T„ 
and Un as jj| 

cos(n6) = Tn(cos{9)) 

and 

sin{n9) = Un-i{,cos{9))sin{9). (19) 



Using these expressions and writing the dot product as a sum of components Ijllll becomes 

3 

5"(p) = Tr,{cos{d))I + iUn-i{cos{e)) ^ Qsm(0)a,. (20) 

The equahties of ifTZJI enable us to rewrite this as 

S^{p) = Tn{cos{P)cos{p' )) I — iUn-i{cos{P)cos{p'))[sin{P)sin{p")ai — sin{j3)cos{p")a2 + cos(/3)sm(p')cr3] (21) 
Using the PauU matrices the matrix polynomial 

S^{p) = 



(22) 



T„{cos{p)cos{p' ) Un-i{cos{[3)cos{p'))sin{(3)cos{p") 

-Un-i{cos{f3)cos{p'))sin{P)cos{p") Tn{cos((3)cos{p')) 

Un^i{cos{j3)cos(p'))cos(j3)sin{p') Un-i(cos{l3)cos{p'))sin{(5)sin{p") 
Un-i{cos{P)cos{p'))sin{f3)sin{p") —Un-i{cos((3)cos(p'))cos{(3)sin{p') 

is obtained. 

The evolution of the quantum walk in momentum space representation given in (|Hl)can also be expressed as 

(/.(n,p)e-"- = 5"(p)</)(0,p). (23) 

lf22|l and © enable this expression to be written as 

</)o(n,p)e-''"" = [Tn{cos{l3)cos{p'))-iUn-i{cos{l3)cos{p'))cos{l3)sin{p')]^^{Q,Q) 
+ [Un~i{cos{(3)cos(p'))sin{l3)cos{p") - iUn^i{cos(f3)cos{p'))sin{l3)sin{p")\^i{Q, 0) (24) 



(/)i(n,p)e~™" = -[C/„_i(cos(/3)cos(p'))sm(/3)cos(p")+iC/„-i(cos(/?)cos(y))sm(/3)sm(p")]^o(0,0) 

+Tr,{cos{(3)cos{p'))+iUn-i{cos{P)cos{p'))cos{l3)sin{p')]^i{Q,0) (25) 

By using the relation 

Tn{x)^Un{x)-xUn-l{x) (26) 

this can be written as 

0o(",p)e""" = [Un{cos{[3)cos{p')) - C/„_i(cos(/3)cos(p'))cos(/3)[cos(p') + ism(p')]]*o(0, 0) 

+ [[Un-i{cos{l3)cos{p'))sin{l3)[cos{p") - szn(p")]]*i(0, 0) (27) 

(/)i(n,p)e"™" = -[C/„_i(cos(/3)cos(p'))sm(/3)[cos(p") +ism(p")]]*o(0,0) 
+ [Un{cos{l3)cos{p')) - C/„_i(cos(/3)cos(p'))cos(/3)[cos(p') +ism(p')]]^i(0,0). (28) 

Inverting the de Moivre formula and moving the global phase term to the right hand side gives the analytic 
expressions 

(/.o(",p) = e'""[C/„(cos(/3)cos(p')) - f/„_i(cos(/3)cos(p'))cos(/3)e*P]*o(0, 0) 

+e'""[C/„_i(cos(/3)cos(p'))sm(/3)e-*P]*i(0,0) (29) 

0i(n,p) = -e™"[C/„_i(cos(/3)cos(y))e^P]*o(0,0) 
+e''^°'[Un{cos{l3)cos{p')) - [/„_i(cos(/?)cos(y))cos(/3)e-'P]^i(0, 0) (30) 

for the general momentum space state functions for a one dimensional quantum walk at time n. 
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Figure 1: Momentum Space Density functions for /3 ; 
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Figure 2: Momentum Space Density functions for /3 = -I 



MOMENTUM SPACE DENSITIES 



The denisity |(/)o(p : i)P for a = 7 = 5 = 0, ^o(0,0) = 1 and ^i(0,0) = is plotted in figures [H H for /? = f , f 
and ^ and for times t = 10, 30, 50, 70. When /3 is fixed the dominant feature of the time series is an increase in 
oscillation frequency with time. This corresponds to the increase in support of the position space densities with time. 
The effect of increasing /? is to trade a decrease in the constant component of the density function for an increase 
in the oscillatory component. This corresponds to a shift in the position space of probability density from the zero 
region of the walk to the outer edges of the walk. 

The sequences shows that the densities converge to a limit as time increases. They also illustrate the fact that 
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Figure 3: Momentum Space Density functions for (3 • 



the momentum space is an attractive representation in which to derive this Umit because the domain of the wave 
functions is constant, p G [— tt, tt]. This is in contrast to the real space where the domain expands with time. 

VI. CONCLUSIONS 

It has been shown that the momentum space dynamic equations for a quantum walk can be derived using a z 
transform of the position space equations for the dynamic walk. An exponential representation of the momentum space 
time evolution operator was derived by using Lie group theory. This enabled the calculation of general momentum 
space wave functions in terms of Chebyshev polynomials. Some simple calculations of the momentum space probability 
densities illustrate the convergence of the momentum wave functions to a limit as time increases. 
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